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Polymer stretching in random smooth flows is investigated within the framework of the 
FENE dumbbell model. The advecting flow is Gaussian and short-correlated in time. 
The stationary probability density function of polymer extension is derived exactly. The 
characteristic time needed for the system to attain the stationary regime is computed 
as a function of the Weissenberg number and the maximum length of polymers. The 
transient relaxation to the stationary regime is predicted to be exceptionally slow in the 
proximity of the coil-stretch transition. 



1. Introduction 

The ability of polymers to considerably change the large-scale statistics of the ad- 
vecting flow has i mportant practical applications, drag reduction being one of the most 
relevant ones (see lGvr fc Bewersdorfjll99.'j) . Polymers affect the dynamics of the advect- 
ing velocity field only if they are highly elongated. Understanding how a single polymer 
chain is stretched by a random flow is thus the first issue to address in the study of 
hydrodynamical properties of polymer solutions. 

At equilibrium, the radial shape of coiled polymers is spherical due to their entropy. 
When placed in a non-homogeneous flow, polymers are deformed and stretched by the 
gradients of the velocity. The product of the longest relaxation time of polymers and the 
characteristic rate of deformation is called the Weissenberg number Wi. For small Wi 
the entropic force prevails and polymers are in the coiled state. When Wi exceeds a 
critical value the molecules become highly elongated and their extension increases sharply. 
This phenomenon is called coil-stretch transition and the critical Weissenberg number is 
known to be approximately one. 

We investigate the statistics of polyme r extension in the FENE (finite extensible non- 
linear elastic) dumbbell model (see e.g. iBird et allll977t for a review). A polymer is 
described as two beads joined by an elastic spring. The elastic force diverges as the 
elongation of polymers attains its maximum value, R m , and this gives a large-extension 
cutoff. Consequently the stationary regime exists however strong the velocity gradients 
are. Since attention is directed only to the dynamics of a single molecule, the feedback 
on the advecting flow is disregarded. To allow analytical progress, the random flow is 
chosen to have the Batchelor-Kraichnan statistics. This means that the considered flow 
is Gaussian, white in time and linear in space. The Batchelor-Kraichnan model is a 
fully solvable model for passive turbulent transport which can provide useful conne ctions 
between theory and real behaviours fsee iFalkovich. Gawedzki fc Vergassolall200lL for a 
review on the applications to scalar and magnetic fields). The results should be intended 
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as a qualitative description of real polymer dynamics. We derive the complete form of the 
stationary probability density function (p.d.f.) of polymer extension and describe how the 
statistics of polymer stretching changes with increasing velocity gradients. Concerning 
the statistics at finite times, we compute the typical time needed for the system to reach 
the steady state and predict how it depends on the maximum length of polymers and the 
stretching by the flow. Our analysis shows that the coil-stretch transition of polymers is 
characterised by an anomalous dynamics in time. 

Th e coil-stretch transition was predicted in 1974 for shear and h yperbolic flows i De^GennesI 
1974) and has been wid ely studied experimentally for such flows llPerkins. Smith fc Chul 



1997llSmith. Babcock fc Chull999tlHur eraJ2002HBabcock et alhoO$i . In contrast, the 



experimental study of polymer dynamics in random flows is a very recent achievement. 
This is due to the difficulty in generating a flow that is random at scales comparable with 
the size of polymers (abo ut 100 urn). This difficu l ty can be overcome thanks to the elastic 
turbulence discovered bv lGroisman fc SteinberJ l)2000(k the flow of a highly elastic poly- 
mer solution at low Reynolds numbers, but large Wi has all the main properties of fully 
developed turbulence. Therefore, in solutions of sufficiently elastic polymers it is possible 
to excite turbulent motion in exceedingly small volumes. Exploitin g elastic turbulence in 
polymer solutions, iGerashchenko. Chevallard fc Steinberd l|2005|) thus investigated the 
stretching and the deformation of a single DNA molecule in a three-dimensional random 
flow. 

Theoretical studies concerning the co il-stretch transition in random flows focused 
mainly on the Hookean dumbbell model llchertkovll2000l iBalkovskv. Fouxon fc LebedevI 
200ol lCelani. Musacchio fc Vincenzil2005l) . This model is suitable only for the coiled state 
(Wi < 1) since the linear force in principle allows infinite extensions and for large Wi 
polymers can become more and more elongated under the action of velocity gradients. 
For Wi > 1 a stationary p.d.f. of the extension no longer exists and this behaviour was 
conjectured to coincide with the coil-stretch transition. To overcome this oversimplifi- 
cation, the maximum length of polymers must be taken i nto account. One possibility is 
to replace the Hookean force by a nonlinear elastic force. IChertkovl Ij200fl) obtained the 
large-value tail of the stationary p.d.f. of the extension for a general anh armonic force; 
such approximate analysis was subsequently applied bv lThiffeaultl l|2003l) to the FENE 
model. Here we exactly derive the complete statistics of polymer stretching within the 
context of the FENE dumbbell model at general Wi. 

Concerning the s tatistics at f i nite ti mes, preliminary results were obtained for the 
Hookean model bv ICelani et al\ l)2005(l . There, t he relaxation time to the stationary 
regime could be defined only in the coiled state. ICelani et all l)2005|) thus derived the 
behaviour of the relaxation time for very small Wi and observed a divergence for Wi = 1; 
this suggested a critical behaviour close to the coil-stretch transition. We present the first 
prediction of the complete dependence of the transient relaxation time on Wi and R m 
with the more realistic FENE model. 

The rest of this paper is organised as follows. In f0we introduce the model and present 
the main results. The stationary p.d.f. of the elongation and the transient relaxation time 
are computed in and in 21 respectively. In ^3 we discuss the relevance of our results 
for experiments and numerical simulations. 



2. Coil stretch transition 

In elastic dumbbell models a polymer is described as two beads connected by a spring. 
The beads represent the ends of the molecule and their separation is a measure of the 
extension. The beads experience: (a) a hydrodynamic drag force modeled by the Stokes 
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law; (b) a Brownian force due to thermal fluctuations of the fluid; (c) an elastic force 
due to the spring connecting one bead to the other. We consider two-dimensional and 
three-dimensional flows indifferently, the dimension of the flow being denoted by d. Since 
in physical applications the elongation of polymers is always smaller than the viscous 
scale of the flow, the dumbbell is assumed to move in a linear velocity field v(r,t) = 
v (t) + r ■ Vt>(<). Inertial effects and hydrodynamic interactions between the beads are 
neglected. Consequently, the separation v ector between the beads, R, evolves according 
to the stochastic differential equation fsee lBird et al\[l97l\l 



dR= (R-Vv)dt- EVQ dt+ ^L dW ^ (2.1) 

T V T 

where Rq is the equilibrium length of the polymer, t is its relaxation time in the absence of 
flow, and W is a d-dimensional Brownian motion which accounts for thermal noise. In the 
FENE dumbbell model, the elastic force F(R) takes the form F(R) = R/(l - R 2 /R 2 n ), 
where R m denotes the maximum extension o f the molecule. In physical applications the 
ratio R m /Ro usually lies between 10 and 100 (|Bird et a/.lll977|) . The length of the vector 
R is a measure of the extension of the polymer. 

Within the Kraichnan model v(x,t) is a statistically ho mogeneous Gaussian field with 
zero mean and second-order correlation l|Kraichnanlll96^ 

(v i (x,t)v j (x + r ) t')) = D l3 (r)S(t-t'). 

In the so-called Batchelor regime the flow is assumed to be smooth in space. If we 
further impose incompressibility and statistical i nvariance with respect to reflections and 
rotations, the tensor Dij(r) must take the form l|Monin fc Yaglomlll975(l 

Dij(r) = D S zj - D x [(d + l)S zj r 2 - 2^-], 

where D Q represents the eddy diffusivity of the flow and Di determines the intensity 
of turbulent fluctuations. In random flows the Weissenberg number can be defined as 
Wi = At, where A is the maximum Lyapunov exponent of the flow, that is the average 
logarithmic growth rate of nearby fluid particle separations. The maximum Lyapunov 
exponent of the Batchelor-Kraichnan flow has asy mptotically a Ga ussian p.d.f. with 
mean value A = D\d(d — 1) and variance A = 2\/d l|Kraichnanlli97l . 

2.1. Stationary regime 

The statistics of polymer elongation is described by the p.d.f. of the norm of R averaged 
over velocity realizationsf : V(R,t) = J (P(R,t))R d ^ 1 dft, where d£l denotes integration 
over angular variables. When the flow v has the Batchelor-Kraichnan statistics 7 , (i?,<) 
obeys a one-dimensional Fokker-Planck equation with nontrivial drift and diffusion co- 
efficients (see Q. Under reflecting boundary conditions (the probability does not flow 
outside the domain of definition) the system reaches a steady state for all Wi\. The 
stationary p.d.f. of the elongation, V s t{R) = lim^oo V(R, t), has the form (see 33 

j».m-*r«»(i + T ^) (1-53-) o< fl <«,„ (2.2) 

where h = [2(Rq/ R 2 n + Wi/d)]^ 1 and N is the normalization coefficient (see l|3.5|) below). 
The stationary p.d.f. is shown in figurenfor different Wi. For small elongations compared 

f Because of the statistical homogeneity of v, the average p.d.f. of the elongation does not 
depend on the point of application of the vector R. 

J This should be contrasted with the Hookean model where the stationary regime does no t 
exist for Wi> 1 fsee lChert,koMl200ntlBalkovskv et aZ.l2000l: IThiffeaultteOOfl: ICelani et alMQ 
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Figure 1. Stationary p.d.f. of polymer elongation for the three-dimensional FENE model at 
different Weissenberg numbers Wi (R m /Ro = 50). 



to the equilibrium length, V s t{R) scales as this result holds for a general elastic 

force since the left tail of V s t comes from the events where the elastic force dominates and 
equation (|2.1I) reduces to a <i-dimensional Langevin equation (for physically meaningful 
elastic interactions F(R) should scale as R for R — > 0). For intermediate extensions, 
i?o *C R <C Rm, the stationary p.d.f. is proport i onal to the power law i^- 1 - 2 ' 1 m 
accordance with the prediction of iBalkovskv et all l)200fj() . For large elongations V s t(R) 
scales as (R m — R 2 ) h and vanishes for R = R m . In practic al applicat i ons. R n/R m <C 1, 
the exponent h is approximatively d/(2 Wi), as predicted bv lThiffeaultJ l)2003(l . Obviously, 
whe n Rm — > oo and Wi < 1, V s t{R) tends to the stationary solution of the Hookean model 
(see ICelani e.t oj]l2005h . 

The maximum of the p.d.f., i?*, determines the fraction of polymers which are highly 
stretched. The graph of i?* as a function of Wi is shown in figure [21 When Wi is smaller 
than one, R+ is of the order of Ro and most of polymers have the coiled equilibrium 
configuration. With increasing Wi the most probable elongation R+ grows slowly until Wi 
exceeds d/(d — 1). Then, a sharp transition occurs to a strongly elongated state. This 
can be appreciated from the behaviour of the first-order derivative of i?* as a function 
of Wi (figure EJ. As Wi becomes very large, i?* approaches R m - The same analysis holds 
for the average extension /j,, apart from the fact that it starts increasing for a smaller 
Wi and its limiting value is |i? m (see figure 0). It is worth noticing that the coil-stretch 
transition becomes sharper and sharper with increasing R m (not plotted). 

The normalized r.m.s. value of the extension, a/^i, a 2 = J(R — n) 2 Vst(R) di?, is rep- 
resented in figure It increases at low Wi until it reaches a maximum value; then a is 
compensated by the sharp increase in \x and at large Wi the rescaled r.m.s. eventually 
relaxes to the constant value l/\/T5. 

The skewness y = [f(R — ^) 3 V s t(R) dR]/a 3 is positive for small Wi and becomes neg- 
ative at large Wi (figure 13 accordingly with the qualitative behaviour of the stationary 
p.d.f. (figure 0. The maximum of skewness in the neighbourhood of the coil-stretch 
transition can be easily understood as follows. At low Wi the p.d.f. is peaked at Rq and 
the skewness is positive. With increasing Wi the right tail starts raising, but /i is still 
of the order of i?o : the skewness, therefore, increases and achieves is maximum value. 
Beyond the coil-stretch transition /j, starts moving towards the maximum extension and 
the skewness decreases until it becomes negative at large Wi, that is when the p.d.f. has 
a long left tail. The limiting value of y for Wi — > oo is — | y/b/3. 
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Figure 2. Left: Most probable rescaled elongation R+fRm and average rescaled extension 
n/Rm as functions of the Weissenberg number Wi (d = 3, R m = 50, Ro = 1). Right: First 
derivative of R*/ ' R m with respect to Wi. 





Figure 3. Left: Normalized root mean square a/fiasa function of the Weissenberg number 
Wi (d = 3, R m = 50, Ro = 1). Right: Skewness y vs Wi for the same values of the parameters. 



2.2. Relaxation to the stationary regime 

We now turn to the time dependence of the p.d.f. of the elongation. Starting from an 
initial condition peaked at i?o, the system relaxes to the stationary regime described 
by H2.2fl . The time needed to reach the stationary regime, T, is solution of a transcendental 
equation which involves continued fractions (see J1J). 

For small Weissenberg numbers, < Wi < d/(d + 4), the transient relax ation time T 
behaves according to the prediction of the linear model ijCelani et aZ.ll2005^ : 

T/t= i[l- Wi(d + 2)/d]~ 1 (2.3) 

independently of R m (see inset in figure 

In the proximity of the coil-stretch transition T displays a maximum as a function of 
the Wi. The relaxation is exceptionally slow in this range of Wi because the stationary 
regime results from the competition between the coiled state and the highly stretched 
state. The position and the value of the maximum relaxation time T max depend on the 
cutoff R m (figure 0}. As the maximum allowed extension of polymers increases, T max 
is closer and closer to Wi — 1 and grows; at large R m the FEN E model should ind eed 
match the Hookean model, where T diverges as Wi tends to one llCelani et oJ.ll2005l) . 

For very large Weissenberg numbers the stretching time is small compared to r and 
the molecules are expected to rapidly reach the highly stretched configuration. Hence, T 
vanishes as Wi tends to infinity. A numerical fit shows that T scales as Wi -1 at large Wi. 

In the next sections we explicitly derive (|2.2|l and the equation for T. 
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Figure 4. Left: Rescaled time of relaxation to the stationary regime, T/t, as a function of 
the Weissenberg number Wi for three different values of Rm/Ro (d — 3). The inset shows 
the linear-model approximation given by equation 12.311 (dashed line) for R m /Ro = 50: the 
agreement is good up to d/(d + 4) ~ 0.4. Right: Dependence of the maximum rescaled relaxation 
time T m ax/T on the maximum relative extension of polymers R m /Ro. 



3. Fokker-Planck equation 

For a fixed realization of the velocity field the p.d.f. of the end-t o-end vector, P(R, t), 
satisfies the Fokker-Planck equation associated with (|2.1() (see e.g. lRiske"nlll989|) : 



d t P + div fl 



R 



= ^V 2 R P (3.1) 



T 



To obtain an equation for V(R, t), we have to average the above equation over the velocity 
realizations and integrate the result over angular variables. The terms of the type (viP) 
in general do not lead to a closed form for the mean p.d.f. and a closed equation cannot 
be deduced from 13.1fl . The Gaussianity and the ^-correlation in time of the Batchclor- 
Kr aichnan model provide an exact closure. Exploiting the Novikov-Furutsu formula (see 
e.g. lKlvatskin. Wovczvnski k, Gurariell996|K we obtain: ((ViVj)P) = —Cijki du e [Rk(P}] 
with Cijkt = Di[(d + l)SikSjg — SijSki — 5u5jk\- We can thus derive from a one- 

dimensional Fokker-Planck equation for P(R,t): 

d s V(R, s) - -d R [A(R)V(R, s)} + d% [B{R)V(R, s)}, (3.2) 

where the time has been rescaled with r, s = t/r, and the drift and diffusion coefficients 
have the form 

A(R) = ( d + ^ WiR - F(R) + (d - 1)% and B(R) = ^R 2 + R%. (3.3) 
d R d 

The coefficients A and B are time independent due to the stationarity of the advecting 
flow. If i?o is set to zero, then equation ]3.2I) r educes to t he approximate equation for 
the large-value tail of the p.d.f. derived by Chcrtkovl (l200(tl . 

To solve Ij3.2|l we impose reflecting boundary conditions, that is the probability current 
associated with the solution, J(R,s) — A(R)T'(R,s) — 3r [B(R)T'(R, s)], vanishes in 
R = and R = R m for all s > 0. This means that there is no flow of probability 
through the boundaries of the domain. Un der these conditions the stationary p.d.f. of 
the elongation takes the form (Riske HTflsl 



q r i-R 

™ = BjR) 



A(x)/B(x)dx 



(3.4) 



where the constant C and the lower integration limit R\ are fixed by the normalization 
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condition. The above formula holds for a general elastic force of the form F(R) = f(R)R. 
Replacing the force of the FENE model into 1)3.4)1 . we thus obtain 1)2.2)1 with 

2T(d/2 + h+l) , 
Rt T(d/2)T(h + 1) 2 F l (d/2, h; d/2 + h + 1; - WiR 2 m /dR 2 ) ' 

The function 2-F1 in the normalization coefficient denotes the hypergeometric function. 



4. Relaxation time 

The time dependent sol ution of the Fokker-Planck 1)3.2)1 can be obtained by separation 
of variables l|Riskenl ll9891. In other words, V(R, s) can be sought in the form 

00 

V(R, s) = V st {R) + c k e-^ s Pk (R) (4.1) 
fe=i 

where the coefficients Ck are fixed by the initial condition V(R,0) and Pk{R) are 
respectively the eigenvalues and the eigenfunctions of the ordinary differential equation 

^_[B(R)p k (R)} - ±[A{R)p k {R)]+^ k p k {R) = 0. (4.2) 

The above equation should be solved with reflecting boundary conditions: Jfe(0) = 
liniR—,^ Jk(R) — 0, Jfc being the probability current associated with the eigenfunc- 
tion p k - It can be shown that the eigenvalues fik are real and non- negative, V s t(R) 
belonging to the eigenvalue ^0 = fe.iskenl H9891. As we will see, the eigenvalues form a 
countable set and may be arranged in ascending order: < fix < fi2 < . . . . The reciprocal 
of /ii , therefore, is the time of relaxation to the stationary regime rescaled by r. 

Equation 1)4. 2JI is a second-order linear differential equation with four regular singu- 
larities in the complex plane. By the change of dependent and independent variable 
z = (R/R m ) 2 , Pk(z) = z^ 1 )/ 2 (1 — z) h Wk(z), such equation can be transformed into a 
standard Heun equation for the function Wk(z) (see lRonveauxlll995l for a review): 

d 2 w k ^ h H 8 | e \ dw k a[3z-q Wk ^ Q (43) 

dz 2 \z z — 1 z — aj dz z(z — l)(z — a) 

where 



dz 2 \z z — 1 z — aj dz z(z — l){z — a) 
d R 2 _ d / [ik 



a = — 



WiR 2 2 V 2Wi 




h+ i+u d ( d . 



4 4\/ V Wi 

7 = d/2 S = h e=l + h 

with h = [2(i?o/i? 2 , l + Wi/d)] -1 . Reflecting boundary conditions for pk map into the 
following limiting conditions for Wk'. 

lim z" / - 1 w k {z) = and lim(l - z) s - 1 w k (z) = 0. (4.4) 

z^O z— >1 

The Heun equation is the general Fuchsian equation with four singularities. In the stan- 
dard form 1)4.3)1 the singular points are 0, 1, a, 00. Let zq be a generic singularity of 1)4.3(1 . 
From the theory of Fuchsian equations, the local behaviour of Wk{z) near zq is specified 
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by th e characteristic exponents pi, p2 associated with zq (see e.g. IWhittaker fc Watsonl 
If pi — p2 is not integer, in a neighbourhood of zq which excludes the nearest other 
singularity Wk(z) can be written in the form b\(z — Zo) Pl ip(z — zq) + bi(z — zo) P2 ip(z — zo), 
where bi, bi are constant and tp, ^ are analytic functions such that ip(zo) ^ 0, ip(zo) =/= 0. 
If pi — p2 is integer and pi > p2, the function i\) can be no longer analytic in zq and 
involve the function log(z — z ). 

The singularity z — has characteristic exponents and 1 — 7; the singularity z = 1 
has characteristic exponents and 1 — S. In physical applications we can exclude the 
situation where 5 is integer. On the contrary, 1 — 7 is zero when d = 2. 

Consider first the case d = 3, where there are not logarithmic singularities in z = 0. To 
fulfill conditions (|4.4(l . Wfc must be simultaneously a local solution about z = and z = 1, 
in both cases belonging to the exponent 0. Such a solution is called a Heun function of 
class I relative to the points and 1, and exists only for a countable set of values of q and 
hence of pk (see lRonveauxlll9 95^ . The condition for the aforeme ntioned Heun function 
to exist leads to a transcendental equation for the eigenvalues pk (|Er delvilj 1 944^) : 

Li — — 



where 

(i + a- l)(i +p- l)(i + 7 - l)(i + w - 1) 



Li = q + ai(i + w) 



(2i + w-l)(2i + w-2) 

ei(i + o;)(7- <5) + + oj) + a^][2i(« + to) + 7(0; - 1) 



Mi 



(2i + w- 1)(2«- 
(i + 1)(« + w - a + l)(i + w - + l)(i + S) 



(2i + uj + l)(2i + LO + 2) 

with to = 7 + S — 1. The rescaled relaxation time T/t is then the reciprocal of the lowest 
non-zero solution of Q4.5|) . In the case d = 2 the conclusions are unchanged since the 
solution involving a logarithm in the neighbourhood of 2 = should be discarded. 

We solved (|4.5() numerically: the continued fraction was computed by the modified 
Lentz's method and the first non- zero solution was evaluated by the root false position 
method (see e.g. lPress et aZ.l[l993|) . 



5. Summary and discussion 

The goal of the paper was to investigate polymer stretching in a turbulent flow within 
the context of a fully solvable model. The statistical features of the Batchelor-Kraichnan 
flow allow us to derive the complete form of the stationary p.d.f. of polymer elongation 
for a general elastic force. When specializing to finitely extensible polymers we recover 
the main properties of polymer dynamics in real turbulent flows and compute the time 
of relaxation to the stationary regime. 

It should be noted that the velocity field we consider is statistically isotropic. Together 
with the ^-correlation in time, this is a key assumption in order to derive a fully analy tical 
solution of the problem. In the experimental setup of iGerashchenko et all l)2005|) the 
elastic turbulent flow is superimposed to a mean shear flow. The long-time statistics 
of polymer extension i n the presence of a mean shear ha s been recently consid e red b y 
Chertkoy etaJ\ <|20r)5l) : ICelani. Puliafito fc TuritsTnl ^OO.^ i lPuliafito fc Turitsvnl l^nOSlk 
"Turitsvul 4200,f) . 
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The main result of our study is the behaviour of the time of relaxation to the steady 
state as a function of Wi. At low Wi the transient relaxation time is an increasing 
function of Wi, it is maximum close to the coil-stretch transition, and eventually tends 
to zero with increasing Wi. Knowing the dependence of the transient time on Wi is 
relevant both for numerical simulations and experiments. For example, in the former 
case, the time required for uncorrelated polymer chains that are suddenly exposed to 
the same flow to correlate is (implicitly) related to the sharpness of stress gradients one 
can expect in the flow. Hence, the prediction of the transient time in our study is useful 
to estimate the required grid spacing to fully resolve those gradients (L. Collins 2005, 
private communication). In the latter case, the fact that the transient relaxation time is 
especially long just below the coil-stretch transition implies that within such range of Wi 
experimental measures are more sensitive to statistical fluctuations. 

Experiments concerning the transient relaxation to the stationary regime can inves- 
tigate the time dependence of the conditional p.d.f. V(R, t\R , 0), which corresponds to 
the initial condition peaked at the equilibrium size: V(R, 0|i?o, 0) = 5(R — Rq). Such 
initial condition can be fixed experimentally as follows (A. Celani 2005, private commu- 
nication). The p.d.f. of the extension is constructed by following the motion of different 
polymer molecules and collecting R(t) for each molecule: one should then start counting 
time only when the length of the corresponding polymer is approximatively Rq- This is 
equivalent to selecting the initial state where all molecules have the equilibrium extension, 
i?(0) = R . 

As for the transient relaxation time, this can be measured directly from the time 
behaviour of the conditional moments of the extension: R n (t) = J R n V(R, t\Ro, 0) dR, 
where n is a positive integer. The conditional p.d.f. can indeed be expanded as in Ij4.1|l 
and the order of series and integral can be interchanged in the definition of R n d ue to the 
integ rability of R n pk(R) and the uniform convergence of series l|4.1|l (see e.g. ISmirnovl 
Il984l) . Therefore, all the moments of the extension converge to their stationary value 
with the same rate as the p.d.f. of the extension. 

To conclude, we believe that the results obtained for the nonlinear dumbbell model 
are relevant for the comprehension of polymer dynamics in turbulent flows at any Weis- 
senberg number. Moreover, we think that our study may stimulate new experiments 
directed to investigate the transient relaxation to the stationary regime. 
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